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ABSTRACT 

We reanalyze the photoevaporation problem of subgalactic objects irradiated by ultraviolet background 
(UVB) radiation in a reionized universe. For the purpose, we perform three-dimensional radiation 
smoothed-particle- hydrodynamics (RSPH) calculations, where the radiative transfer is solved by a direct 
method and also the non-equilibrium chemistry of primordial gas including H2 molecules is incorporated. 
Attention is concentrated on radiative transfer effects against UVB for the formation of subgalactic 
objects with T^-„ < lO'^K. We consider the reionization model with Zreion ~ 7 and also the earlier 
reionization model (zrcion ~ 17) inferred by the WMAP. We find that the star formation is suppressed 
appreciably by UVB, but baryons at high-density peaks are self-shielded even during the reionization, 
forming some amount of stars eventually. In that sense, the photoevaporation for subgalactic systems 
is not so perfect as argued by one-dimensional spherical calculations. The final stellar fraction depends 
on the collapse epoch and the mass of system, but almost regardless of the reionization epoch. For 
instance, a few tenths of formed stars are born after the cosmic reionization in Zrcion ~ 7 model, while 
more than 90% stars are born after the reionization in the WMAP reionization model. Thus, effects of 
UVB feedback on the substructure problem with a cold dark matter (CDM) scenario should be evaluated 
with careful treatment of the radiative transfer. 

The star clusters formed at high-density peaks coalesce with each other in a dissipationless fashion in 
a dark matter potential, resultantly forming a spheroidal system. As a result, these low-mass galaxies 
have large mass-to-light ratios such as observed in dwarf spheroidals (dSph's) in the Local Group. 



Subject headings: galaxies: formation 
hydrodynamics 



galaxies: dwarf — radiative transfer — molecular processes 
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1. INTRODUCTION 

In the context of cold dark matter (CDM) cosmology, 
the first generation of objects should have the mass of 
~ IO^Mq and form at redshifts of 20 < z < 50 (Tegmark 
et al. 1997; Fuller & Couchman 2000). At later epochs, the 
first objects are assembled into larger systems in a hierar- 
chical fashion to form dwarf or normal galaxies. On the 
other hand, it is thought that the universe was reionized 
at z > 5, because the absorption by the neutral inter- 
galactic medium (IGM) is observed to be feeble at red- 
shifts z < A. The reionization redshift is inferred to be 
2icion ~ 6 from the spectra of high-z quasars (Djorgovski et 
al. 2001, and references therein), 6 < Zroion ^ 10 from sim- 
ulations on Lya absorption lines (Umemura, Nakamoto, & 
Susa 2001), or 11 < Zroion ^ 30 from the recent results of 
WMAP (Kogut et al. 2003). Thus, it is anticipated that 
the UVB directly influences the assembly phase of galax- 
ies. The UVB has several important physical impacts on 
the galaxy formation. One of the most important effects 
is the photoheating. If a gas cloud is irradiated by UVB, 
the gas temperature is raised up to T ^ 10"* — lO^K, so 
that gaseous systems with virial temperature of T < 10^ 
K can be evaporated owing to the enhanced thermal pres- 
sure (Umemura & Ikeuchi 1984; Dekel & Rees 1987; Efs- 



tathiou 1992; Babul & Rees 1992; Thoul & Weinberg 1996; 
Ferrara & Tolstoy 2000; Gncdin 2000c; Shapiro & Raga 
2001; Shaviv & Dekel 2003). Such an effect by UVB 
may reconcile the paradox that low mass galaxies are over- 
produced in CDM cosmology, compared with observations 
(White & Frenk 1991; Kauffmann, White, & Guiderdoni 
1993; Cole et al. 1994; Moore et al. 1999; Binney, Ger- 
hard, & Silk 2001). The criterion for photoevaporation has 
been derived from detailed hydrodynamic calculations by 
Umemura & Ikeuchi (1984) and Thoul & Weinberg (1996). 
However, if the systems are self-shielded against the UVB, 
the criterion for photoevaporation is completely changed. 
The self-shielding comes from the radiative transfer effects 
of ionizing photons. Tajiri & Umemura (1998) have de- 
rived the self-shielding criterion by solving full radiative 
transfer for spherical clouds, to find the critical optical 
depth of 2.4. Also, if the system is self-shielded from UVB, 
the gas can cool below lO'' K by H2 cooling (Susa & Ki- 
tayama 2000). The complete radiation- hydrodynamic cal- 
culations have been done for spherical clouds by Kitayama 
et al. (2001), where hydrodynamics, radiative transfer, 
and primordial gas chemistry including H2 molecules are 
self-consistently incorporated, and thereby the criteria for 
cloud collapse have been derived depending on the red- 
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shifts. Very recently, pioneering approaches with three- 
dimensional hydrodynamics have started on this isuue in 
the context of CDM cosmology. They are, for example, 
the simulations by Ricotti, Gnedin, & Shull (2002a, b) 
with radiative transfer based on the optically thin vari- 
able Eddington tensor (OTVET) approximation, and by 
Tassis et al. (2003) with optically-thin approximation. 

If the nonlinear evolution of density fluctuations pro- 
ceed in a hierarchical fashion, two major effects by radia- 
tive transfer arc expected. One is the direct effect, that 
is, the self-shielding of density peaks that collapse prior to 
the reionization (Gnedin 2000c; Nakamoto, Umemura, & 
Susa 2001; Razoumov et al. 2002; Ricotti, Gnedin, & 
Shull 2002a, b). The other is the enhancement of H2 for- 
mation in relatively massive fluctuations which are once 
ionized and then self-shielded in the course of mass accre- 
tion. The second eff'ect has been pointed out by many au- 
thors so far (Kang & Shapiro 1992; Corbelli, Galli & Palla 
1997; Susa & Umemura 2000). These two effects play a 
substantial role for the galaxy formation under UVB, es- 
pecially if the star formation process is taken into account 
in collapsing clouds, because star formation activities in 
galaxies appear to be strongly correlated to cold HI gas 
with a few xlO^K (Young & Lo 1997a,b). Therefore, the 
effects of radiation transfer of ionizing photons should be 
carefully treated to elucidate the galaxy formation process 
under UVB. In this paper, we perform 3D SPH calcula- 
tions, incorporating radiation transfer and primordial gas 
chemistry. Here, we focus on relatively low mass galax- 
ies, where the dynamics is likely to be strongly subject to 
UVB. The evolution of subgalactic systems is significant 
also for the formation of massive galaxies, since the self- 
shielding against UVB can be responsible to overall SFH 
in the galaxies and therefore resultant galactic morphology 
(Susa & Umemura 2000). 

The goal of this paper is to elucidate how the photoe- 
vaporation of subgalactic systems is caused by UVB and 
how the star formation history (SFH) is influenced con- 
sequently. The paper is organized as follows. In §2, the 
numerical methods are described. In §3, the results of 
simulations are presented. In §4, the SFH is discussed in 
relation to the inferred SFH in dwarf galaxies. §5 is de- 
voted to the conclusions. 

2. NUMERICAL METHODS 

We have performed numerical simulations on the single 
halo collapse with SPH particles (2^^) and dark matter 
particles (2^^). The halo mass ranges from Afhaio ~ 6 x 10^ 
to 6 X 10^ Mq. The initial conditions are set by COSMICS. 
The radiative transfer of the UV radiation fleld is solved 
assuming a single source outside the simulation box, cou- 
pled with the detailed chemistry of primordial gas. We 
also have taken into account the numerical "star forma- 
tion" . Following subsections are devoted to describe some 
details of the scheme. 

2.1. Hydrodynamics and Gravity 

Hydrodynamics is calculated by Smoothed Particle Hy- 
drodynamics (SPH) method. We use the version of SPH 
by Umemura (1993) with the modification according to 

^ In this case, the fraction of electrons can be an order of magnitude s 
by metals becomes an order of magnitude smaller, since the cooling is i 



Steinmetz & Miiller (1993), and also we adopt the particle 
resizing formalism by Thacker et al. (2000). The gravi- 
tional force is calculated by a special purpose processor for 
gravity calculation, GRAPE (Sugimoto et al. 1990), which 
can also accelerate an SPH scheme (Umemura ct al. 1993). 
Here, we used the latest version, GRAPE-6 (Makino 2002), 
which has the speed of ITflops. Eight GRAPE-6 boards 
are combined with CP-PACS, which is a massively par- 
allel supercomputer in University of Tsukuba, composed 
of 2,048 processing units, with theoretical peak speed of 
614 Gflops (Iwasaki 1988). This system is called the Het- 
erogeneous Multi-Computer System (HMCS) (Boku et al. 
2002). The hydrodynamics, radiative transfer, and chem- 
ical reactions are solved with CP-PACS. The originally 
developed software for HMCS can be easily applied to any 
massively parallel computer connected through TCP-IP 
with GRAPE-6 boards. Actually, we also used an alpha- 
cluster and a PC-cluster instead of CP-PACS as the host. 

The softening length for gravity is taken as 60pc for all 
SPH and CDM particles. Also, the minimal size of SPH 
particle is introduced so as to prevent the simulation from 
stopping owing to very short local timescales. We have 
tested this scheme by the standard Sod's shock tube and 
Evrard's collapse problems, and the results well reproduce 
the known solutions. 



2.2. Thermal Processes 

The non-equilibrium chemistry and radiative cooling for 
primordial gas are calculated by the code developed by 
Susa & Kitayama (2000), in which H2 cooling and re- 
action rates are taken from Galli & Palla (1998). The 
H2 cooling rate induced by the collision with H atoms is 
plotted in Figure 1 for high (10^cm~^) and low density 
(0.1 cm~^) cases. The cooling functions by HoUenbach & 
McKee (1979) and Lepp & Shull (1983) are also plotted 
for reference. They agree with each other at the high den- 
sity limit (i.e. LTE), whereas they disagree significantly 
with each other at the low density limit. We also assessed 
the contribution of metallic cooling, and compared it to 
H2 and Lya cooling in Figure 2 for the temperature below 
lO^K. The cooling rate by metals is evaluated by the fitting 
formula in Dalgarno & McCray (1972), assuming that the 
relative abundance between the heavy elements is the same 
as the solar neighbourhood. The fractions of electrons and 
H2 are assumed to be the values behind a shock with the 
velocity of Vs > 30kms~^, since the initially ionized gas 
traces a similar path as the gas behind such a shock (Kang 
& Shapiro 1992; CorbelU, GaUi & Palla 1997; Susa & 
Umemura 2000). We also plot H2 cooling function assmn- 
ing the H2 fraction of 10^^, which is favored for the clouds 
with virial temperature of lO^K < Tvir < 10''K(Nishi & 
Susa 1999). 3 As shown in Fig.2, for T < 8 x 10^ K, the 
dominant cooling mechanism is dependent on the metal- 
licity. 

In nearby dwarf galaxies, the metal cooling is likely to 
be important, because the observed metallicity is Z/Zq ~ 
10~^ — lO"'^. On the other hand, the observations of Lya 
absorption systems at high redshifts imply that the metal- 
licity of high-2 intergalactic medium (IGM) is at a level of 

lUer than the high velocity postshock cases. Hence, the cooling rate 
uced by the collision between electrons and heavy elements. 
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ZjZQ - 10-2 _ 10-3 (Songaila 2001). Hence, the cool- 
ing rate by metals is rather smaller than that by H2, as 
long as lO^K < T < lO^K. In the present paper, we fo- 
cus on the primary star-forming phase in dwarf galaxies 
from metal poor gas with the mctallicity seen in IGM or 
more primordial gas. If one pursues the subsequent recy- 
cling phase of interstellar medium with SN explosions, it 
is definitely requisite to include the metal cooling. Such 
a recycling effect is important for the chemical evolution 
of galaxies, but the issue is beyond the scope of present 
paper and left in future study. 

It should be also noted that the effects by metals can 
play a significant role even at a level of Z/Zq ^ for 
the runaway collapse onto a primordial protostar (Omukai 
2000; Bromm et al. 2001; Schneider et al. 2002), al- 
though here we do not pursue the runaway to a star in 
such dense, cold regions. Based upon these arguments, 
wc neglect the cooling by metals in this paper, and the 
thermal processes are calculated with the chemistry for 
primordial abundance gas. 

2.3. Radiative Transfer 

The photoionization rate and the photoheating rate are 
respectively given by 



hv 
hu 



(1) 



ffuihu — huL)dQ.dv. (2) 



Here nui represents the number density of neutral hydro- 
gen, au is the photoionization cross section, z^l is the fre- 
quency at the Lyman hmit, and Vl is the sohd angle. Iv{s) 
is the intensity of the UV radiation with s denoting the 
distance from the source of ionizing photons. Iv{s) is de- 
termined by solving radiative transfer equation: 



dl. , 

ds 



(3) 



where the reemission term is not included, because we em- 
ploy on-the-spot approximation. The transfer equation 
(3) is solved based upon the method developed by Kessel- 
Deynet & Burkert (2000) , which utilizes the neighbor lists 
of SPH particles to assess the optical depth from a certain 
source to an SPH particle. Here, we assume a single point 
source located very far from the simulated region, and con- 
trol the UV intensity by specifying the incident flux to the 
simulation box as described in section 2.5. It is noted that 
the irradiation from one side can overestimate the effects 
of shadowing, and also internal sources, if they form, might 
play a significant role (Ricotti, Gnedin, & ShuU 2002a, b). 
On the other hand, in Kessel-Deynet & Burkert's method, 
the numerical diffusion of radiation is not negligible. This 
effect tends to reduce the shadowing effect. These points 
should be improved in the future. 

We parallelize the scheme so that wc can implement it 
on massively parallel computer systems such as CP-PACS. 
For the purpose, particles should be classified into the sub- 
sets which are assigned to PUs. The domain decomposi- 
tion is performed as shown in Figure 3. The first step 
is to project the particle positions onto the plane which 
is perpendicular to the light rays of the external single 
source. Since the position of source is assumed to be far 



away from the simulated region, we can approximate the 
light rays to be parallel. The second step is to sort the pro- 
jected positions of particles on the x — y plane. Then, with 
the results of coordinate sorting, we generate subsets with 
an equal number of particles, which arc respectively as- 
signed to computational domains, i.e. PUs. In Appendix, 
the test calculation using this scheme is presented for the 
propagation of ionization front in a uniform media. We 
remark that the present parallelization scheme is not valid 
for the multiple sources, because the domain decomposi- 
tion is based upon the projection of all particles onto the 
plane perpendicular to the source direction. Hence, it is 
necessary to invoke a new technique in order to include 
internal sources in the future analysis. 

Combining equations (1),(2) and (3), we can rewrite the 
ionization rate and the photoheating rate as follows: 

ds J^, J 
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When we need to assess the rates at each grid point (or 
each SPH particle), wc use the following formula: 

fcion {l + 1/2) = i ($i(i) - $i(z + 1)) , (6) 
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Here z -I- 1/2 denotes the grid point at which the rates 
are assessed, and i denotes the boundary between z — 1/2 
and i + 1/2 grids. In the present simulation, these grids 
are generated by the method in Kessel-Deynet & Burkert 
(2000). The above formula has an important advantage 
that the propagation of ionization front is properly traced 
even for a large grid size with optical depth greater than 
unity (r > 1) (Kessel-Deynet & Burkert 2000; Abel, Nor- 
man & Madau 1999). 

2.4. Star Formation Algorithm 

Here the "star formation" algorithm used in this sim- 
ulation is described. At each time step, we pick up the 
SPH particles that satisfy following four conditions: (1) 
div v < 0, (2) p/p > 200, (3) T < 5 X lO^K, and (4) 
> 5 X IQ-^. These conditions look similar to those 
adopted by Cen & Ostriker (1993), but the conditions (3) 
and (4) are different. The first condition guarantees that 
gas surrounding the particle is infalling. In addition, if the 
region is virialized, the density around the particle should 
satisfy the second condition. The conditions (3) and (4) 
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cannot be satisfied unless the region is self-shielded against 
UVB and thereby H2 cooling is effective. These two con- 
ditions are essential for the star formation, although they 
are not hitherto taken into account in the previous simu- 
lations. Also, the conditions (3) and (4) match the obser- 
vations on the HI contents in nearby dwarf galaxies, which 
indicate that the presence of cold (T ~ lO^K) HI gas is 
a good indicator of star formation activities (Young & Lo 
1997a,b). 

The next step is to define the timcscalc of "star forma- 
tion". In order to define the conversion timescale from 
gas to stars, we use the following simple and often used 
formula: 

dp _ dpstar ^^^^ 



dt 

dpstm 

dt 



dt 



(12) 



where tg is the local free-fall timescale, pstai is the mass 
density of star particles, and c* denotes the parameter 
of star formation timescale. Following the above expres- 
sion, an SPH particle that satisfies the conditions for 
star formation is converted to a collisionless particle, after 
At = tg/c^- Here, we use c* = 0.1 for the fiducial model, 
and also investigate = 1 and = 0.01 for scvc;ral c;ases. 
In the previous numerical simulations on the formation of 
disk galaxies, c* is sometimes set to be 0.05 — 0.1 with ap- 
plying the cooling criterion of lO^K, in order to regulate 
both the star formation efficiency and the star formation 
rate (SFR) to match the Kennicutt's low (e.g. Kennicutt 
1998; Koda, Sofue & Wada 2000). However, in the present 
simulation, the star formation efficiency is physically reg- 
ulated in terms of the temperature criterion (3). Hence, 
in this paper, c* = 1 does not mean that all the gas cooled 
to lO^K is converted into stars. Here, c* just controls the 
star formation timescale or SFR. 

2.5. Setup 

The initial particle distributions in a CDM-dominated 
universe are generated by a public domain code, GRAFIC, 
which is a part of COSMICS. Throughout this paper, 
we assume a A-dominatcd fiat universe, with Qm = 0.3, 
Oa = 0.7, Qb/i^ = 0.02, and h = 0.7. Total mass of the 
simulated region is 10® — 10^ Mq, in which a halo collapses 
at 1 < Zc ^ 10. We first generate the positions and ve- 
locities of particles in a cube with the constraint that the 
peak of an overdense region is located at the center of the 
cube. The peak height is controlled so that the overdense 
region collapses at a given epoch. Then, we hollow out a 
spherical region from the cube, so that the sphere contains 
the overdense region and the radius of the sphere agrees 
with the smoothing scale of the overdensity. Then, we at- 
tach a rigidly rotating velocity field to the spherical region 
which corresponds to spin parameter A = 0.05 (Heavens & 
Peacock 1988; Steinmetz & Bartelmann 1995). In each 
run, we define the halo mass Mhaio(z) as the mass within 
the radius Ryir{z). Virial radius Ryii{z) is defined by the 
condition, 

Pdm{< Ryir{z)) = 187rVo(l + zf, (13) 
where pdm{< R) is the averaged dark matter density 
within the radius R measured from the center of mass of 
the whole system, po denotes the averaged density of the 



present universe. We also define the collapse epoch Zc as 
the epoch by which a half of the dark halo mass collapsed 
within the radius Rvii{z). The typical halo mass just after 
the virialization is approximately 60% of the total mass of 
the simulated region. These are tabulated in Table 1. 

For the largest simulations, 2^'^{= 131072) SPH particles 
-|- 2^^" dark matter particles are employed. Also, the case 
study is done for all parameters with a smaller number of 
particles, 2^^ SPH particles -I- 2^^ dark matter particles. 
Finally, we give the Hubble expansion velocity to all parti- 
cles in addition to the already generated peculiar velocity 
fields. 

We also fix the initial abundance of chemical species. 
We have taken the values of the canonical model in 
GaUi & Palla (1998). The evolution of the ultravio- 
let background (UVB) is modeled as follows. The spec- 
trum shape of UVB is assumed to he 1^ oc u~^, and 
the intensity I21, which is the UVB intensity normal- 
ized by 10~^^erg s~^cm~^str^^Hz^^, is assumed to be 
/21 = [(1 -I- z) /Sf hi z <2 and hi = 1 ior 2 < z < 4. 
This dependence is consistent with the UV intensity in the 
present epoch (Maloney 1993; Dove & ShuU 1994) and 
the value inferred from the QSO proximity effects at high 
redshifts (Bajtlik, Duncan, & Ostriker 1988; Giallongo 
et al. 1996), although there is an uncertainty of lO^*^'^ at 
2 < z < 4. As for higher redshift epochs, we employ two 
regimes. One is a regime with Zroion ~ 7, which is provided 
by I21 = cxp [3 (4 — z)] for z > 4 (Umemura, Nakamoto, 
& Susa 2001). Such UVB evolution is suggested by the 
comparison between Lya continuum depression in high-z 
QSO spectra and the simulations of QSO absorption lines 
based on 6D radiation transfer calculations on the reioniza- 
tion (Nakamoto, Umemura, & Susa 2001). In this regime, 
the reionization proceeds from z 9 and is completed at 
z Rs 7. Most of the runs are performed in this regime. 
B\it, for comparison, we also investigate a high-z reioniza- 
tion regime (Zreion « 17) such as recently inferred by the 
WMAP (Kogut et al. 2003). In this regime, the intensity 
of UVB is modeled by I21 = max(exp [3 (4 — z)] , 0.01) for 
4 < z < 17 and I21 — exp [3 (17 — z)] for z > 17, based 
on Nakamoto, Umemura, & Susa (2001). We calculate 
typical two runs in the high-z reionization model. 

3. RESULTS 

In Table 1, the model parameters studied here and the 
basic results are summarized, where Zc is the collapse red- 
shift, Mhaio is the halo mass, T^ir is the virial temperature, 
and (7iD is the line-of-sight velocity dispersion of stellar 
component, ctid(DM) is that of dark matter component. 
These values should be defined as functions of redshift, 
but the tablated values are assessed just after the virial- 
ization. i?haif is the effective radius of the formed galaxy, 
/star is the final stellar fraction, and Mh^\o/L is the mass- 
to-luminosity ratio in solar units. We show the detailed 
results in the following. 

3.1. Reionization Feedback 

Figures 4 and 5 show two typical results, where c, = 0.1 
for both cases. Figure 4 shows the low-Zc (zc — 1.7) and 
less massive case (Mhaio — 6 x IO^Mq). Dots in the fig- 
ure denote the location of SPH particles and numerically 
formed stars. The colors of particles represent the gas 
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temperature (color legend is shown at the bottom) , except 
that blue particles denote stars. At z ~ 66 (1st panel), the 
initial condition is set in a spherical region as described in 
section 2.5. Thereafter, the system expands by the cosmic 
expansion, and also baryon perturbations arc induced by 
imprinted dark matter fluctuations (2nd panel). At z ~ 9, 
first objects with the mass of ~ 10^~'^Mq form at den- 
sity peaks (2nd and 3rd panels). This is consistent with 
the previous results by Tegmark et al. (1997) and Fuller 
& Couchman (2000). At z « 7, the system is rcionizcd 
overall and reheated by UVB (4th panel). Then, a por- 
tion of the cooled gas that is once settled into the first 
objects is evaporated as a pressure-driven thermal wind at 
z w 5 (5th panel). But a significant amount of baryons 
are self-shielded against UVB at high- density peaks even 
during the reionization, so that the photoheating is pre- 
cluded deep inside the peaks and eventually stars form 
there. As a result, small star clusters arc left after the 
reionization epoch (6th panel). The stellar clusters coa- 
lesce with each other in a dissipationless fashion, accord- 
ing as dark matter fluctuations grow hierarchically (7th 
panel). Finally, a subgalactic spheroidal system with a 
high total mass-to-stellar mass ratio (Mhaio/-^star — 20) 
forms. If we assume the stellar mass to luminosity ratio 
of 3 in solar units that comes from a Salpter-like initial 
mass function, then the mass-to-light ratio of this system 
is assessed to be Mhaio/i = 63. 

In the above case, baryons are cooled to form stars only 
at density peaks, and gas in other regions is photoevapo- 
rated, resulting in a thermal wind. But, the evolution is 
quite different for the high-Zc and massive halo model. Fig- 
ure 5 represents the Zc ~ 7.6 and Mhaio — 6 x IQ^Mq case. 
In this case, most of baryons collapse before the reioniza- 
tion. Resultantly, the photoevaporation is less effective 
than the low-^c case, and thus the whole system does not 
lose a significant portion of baryonic matter. Finally, a 
formed spheroidal system has Mhaio/A^star = 7.7, which 
is close to the assumed initial the total mass-to-baryonic 
mass ratio, fl/ilB = 7.35. The mass-to-light ratio of this 
system is assessed to be M^^aio/L = 26. 

3.2. Star Formation History 

In Figure 6, the star formation histories (SFH) are 
shown for typical four runs, which are the low mass 
(6 X 10 ''Mq) case with Zc = 1.7 and 8.1 and the high 
mass (6 x IO^Mq) case with Zc — 1.1 and 7.6. As com- 
monly seen, the SFR is peaked before lO^yr {z > 5) and 
the star formation activity is substantially suppressed by 
UVB after the reionization, z < 5. This indicates that the 
star formation mainly occurs in high-density regions which 
collapsed prior to the reionization. The peak SFR is lower 
in the systems with lower collapsing epochs. However, the 
difference is not so drastic. The peak SFR in » 1 cases 
is roughly a half of that in Zc « 8 cases. If based on the 
criterion for spherical clouds (e.g. Kitayama et al. 2001), 
the clouds should be completely photoevaporated in the 
two cases of Zc ~ 1. But, the present results show that 
high-density peaks which collapse before the reionization 
contribute significantly to the star formation even in such 
low-Zc objects. This is a consequence of the hierarchical 
growth of density fluctuations that is caused by the CDM 
spectrum. 



As another important result, it should be noted that the 
star formation activity in the case with Mhaio — 6 x IO^Mq 
for high-Zc continues even after the reionization. This is 
caused by the depth of gravitational potential. A mas- 
sive galaxy with high-Zc has a deeper gravitational poten- 
tial. Thus, the dark halo can prevent the photoheated gas 
from blowing out into the intergalactic space. Instead, the 
heated gas slowly accretes in the potential of dark halo. 
Eventually, the gas is self-shielded from UVB and cools by 
II2, so that stars form even in UVB. 

In Figure 7, the final fraction of stars, /star, is plotted 
against collapse epoch Zc, where /star is defined by the total 
mass of star particles divided by the total baryonic mass 
at z = 1 (/star = A'/starA-fbaryon)- The Starred pentagons 
represent the results of Mhaio — 6 x lO^M© runs, and the 
pentagons denote those of Mhaio — 6 x 10 ''M© runs, where 
c* is assumed to be 0.1. The stellar fraction is smaller for 
Mhaio — 6 X 1O''M0 cases than Mhaio — 6 x IO^Mq cases. 
This comes from the SFH shown above. The massive ob- 
jects tend to retain more gas components than less massive 
objects in the presence of UVB. 

Additionally, /star decreases monotonously with de- 
creasing Zc for both Mhaio — 6 x IO^Mq and Mhaio — 
6 x 10'' M0 cases. It is found that the dependence of /star 
on 1 -|- is almost linear, and its relation is given by 

/star ^ a(M) (1 + Zc) + b{M) (14) 

where a(6 x 10*^ Mq) = 7.7 x 10"^ a(6 x lOUfo) = 6.9 x 
10-2 , 6(6x10SMq) = 0.26 and 6(6x IO^Mq) = 7.7x10-2. 
We emphasize that even for galaxies assembled after the 
reionization, some fraction of baryons are cooled and form 
stars at high-density peaks. The vertical two lines denote 
the critical redshifts obtained by Kitayama ct al. (2001), 
after which protogalactic clouds are photoevaporated. The 
condition for Mhaio 6 x W'^Mq is denoted by a solid line 
and that for Mhaio =i 6 x IO^Mq by a dotted Une. This 
result shows again that the inhomogeneity in protogalactic 
clouds is crucial for the star formation in them, because 
a portion of baryons collapse and cool before the reion- 
ization, and also some fraction of gas is self-shielded from 
UVB even after the reionization if the gas accretes in dark 
matter potential. 

In Figure 7, we plotted also the results with = 1 
(high SFR runs) and c* = 0.01 (low SFR runs). As seen 
clearly, the results are basically the same as the case of 
c* = 0.1. This is physically imderstood as follows: In 
the present simulations, stars form from overdcnsc regions 
with p > 10-22g cm^^. Then, the star formation timescale 
assessed by the present star formation algorithm (12) is 
given by tsF = tg/c^ < lO'' c-^yr. On the other hand, the 
UV feedback is quite effective at z < 5 which corresponds 
to tn > lO^yr. Thus, the main episode of the initial star 
formation ends before the reionization as far as c* > 0.01. 
Moreover, the density of star forming regions is higher if 
c* is smaller, because further radiative cooling works be- 
fore the gas is converted to stars. Therefore, the local 
free-fall time of star forming region is shorter for smaller 
c* . Consequently, c* is not so important to determine the 
final stellar fraction. However, if we adopt = 0.01, the 
star formation history is modified as a matter of course. 
The epoch of star formation is shifted to later epoch for 
c* = 0.01, as shown in Fig. 8. Especially, in the cases with 
high mass and high-^c, the star formation activity after 
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the reionization is enhanced compared to the c, = 0.1 
case. The relation of this result to the continuous star 
formation history inferred in nearby dwarf galaxies (e.g. 
Mateo 1998) will be discussed later. 

3.3. Radiative Transfer Effect 

Here, we show more specifically the effects of radiative 
transfer in the present simulations. For comparison, we 
perform the control runs in which gas is assumed to be 
optically thin against UV radiation. In Figure 9, the time 
variations of stellar fractions are compared between the 
runs with and without radiative transfer. The parameters 
arc the same as those in Figure 6. For all runs, it is clear 
that the star formation activities are supressed at earlier 
epochs and thoroughly terminated after the reionization in 
optically-thin simulations. In Figure 10, the final stellar 
fractions (/star) are shown for the runs with and without 
radiative transfer. The optically-thin approximation un- 
derestimates the stellar fraction by a factor of 1.5-2 for all 
runs in the z^cion — 7 models, regardless of the mass and 
collapse epoch. The radiative transfer effect is even more 
crucial in Zroion — 17 models as shown in §3.5 below. The 
reduction of /star by dismissing the radiative transfer is 
more than an order of magnitude if ^reion — 17. Hence, 
it is clear that the self-shielding does work to form stars 
effectively during the reionization. 

Also, in the bottom-right panel of Figure 9, where 
-^haio — 6 X 1O*M0 and Zc = 7.6, /star continues to in- 
crease if the radiative transfer is incorporated. Since the 
gravitational potential is as deep as to sustain the pho- 
toheated gas in this case, some gas is self-shielded in the 
course of accretion onto local density peaks, so that stars 
continue to form. In the optically-thin simulation, the 
accreted gas cannot be cooled below lO'^K owing to the 
absence of shielding effect, so that no stars form after the 
reionization. 

3.4. Kinematic Properties 

In Table 1, ID velocity dispersion of the formed stars 
(ctid) and effective size i?haif of the resultant spheroidal 
system arc listed for various runs. Here, the effective size, 
-Rhaif ) is defined so that a half of stars are contained within 
the radius. In the low-^c case shown in Figure 4, gas can 
cool at density peaks before the reionization and form com- 
pact star clusters. The star clusters coalesce in a dissipa- 
tionless fashion to form a single spheroidal system eventu- 
ally. On the other hand, in the high-^;c case as shown in 
Figure 5, the formation process of the final stellar system 
is rather different, because the star formation proceeds not 
only at the density peaks but also after the collapse of the 
whole system, which takes place after the reionization. In 
order to clarify the difference quantitatively, we plot the 
time variation of the ratio of the total kinetic energy by 
random motion (-Bran) to the total rotation energy (E^ot)- 
Here, Ej-ot is defined as the summation of the rotational 
energy of baryonic particles (i.e. SPH particles and star 
particles) which is contained in a sphere with radius of 
Ryir{z) from the center of gravity; 

-E^rot = X! 2™''^™'''' ^^^^ 

n<TU-„{z) 

where i^rot.i is the angular velocity of *-th particle, i^ran is 
defined by the rest of the kinetic energy, that is, i^ran = 



ii^kin — E,-ot, where i?kin is the total kinetic energy. When 
the velocity field is isotropic, the ratio E^ot / E^an should be 
1/5 according to the present definition. Figure 11 shows 
the time evolution of the ratio -Erot/^'ran for typical four 
runs after the collapse epoch. As shown in this figure, 
E'rotZ-E'ran is wcU bclow Unity for all the cases after Zc- 
This means that these systems are not rotationally sup- 
ported. However, if we take a closer look, the ratio for the 
high-mass and high-^c case (Mhaio — 6 x lO^M©, Zc ~ 7.6) 
is somewhat larger than the others. This reflects the fact 
that a part of gas forms a rotating disk before it is con- 
verted into stars. On the other hand, for the other cases, 
gas is lost after the reionization era {z <7), and no rotat- 
ing disc forms. 

3.5. Effects of Early Reionization 

Here, we attempt to evaluate the effects of early reion- 
ization as inferred by the results of WMAP (Kogut et al. 
2003). We perform two runs in the high-z reionization 
regime (zreion ~ 17) described at the end of §2.5. The 
both runs are high-Zc models, where Mhaio — 6 x IO^Mq 
with Zc = 8.1 or Afhaio =^ 6 x lO^M© with Zc = 7.6. The 
star formation history and the time evolution of stellar 
fractions are shown in Figure 12. First, for the former run 
(A^haio — 6 X lO'^'AfQ), the effects of early reionization feed- 
back are noticeable, and the fraction of stellar component 
becomes roughly four times smaller than the fiducial reion- 
ization regime (zreion ~ 7). Nonetheless, some fraction of 
baryons are cooled to form stars even at z < 17, owing 
to the effects of self-shielding. This is partly because the 
intergalactic medium is not perfectly transparent to UVB, 
even if the UV intensity is strong enough to reionize the 
universe (Nakamoto, Umemura, & Susa 2001). On the 
other hand, for the latter run (Mhaio — 6 x lO^M©), the 
final stellar fraction is smaller than the fiducial model, but 
the effect is not so drastic as the former case. In this case, 
the final fraction of stars is ~ 20% smaller than the fiducial 
model. 

We also stress that the effect of radiation transfer is very 
important for such early reionization model. In the lower 
panels of Figure 12, the results from the optically-thin sim- 
ulations are also plotted by dashed curves. It is quite clear 
that the star formation is almost prohibited if the optically 
thin is assumed, where the final stellar fraction is reduced 
by an order of magnitude. The results of more systematic 
investigation will be reported in a forthcoming paper. 

4. DISCUSSION 

4.1. Substructure Problem 

Recent numerical simulations predict that subgalac- 
tic dark halos are overabundant compared to dwarf 
spheroidals observed around our Galaxy (Moore et al. 
1999). So far, the possibility that this substructure prob- 
lem might be reconciled if the star formation is signifi- 
cantly suppressed by UVB in subgalactic halos has been 
discussed extensively. However, the present simulations 
have shown that the suppression by UVB is not complete. 
Especially, even in the early reionization regime, the pho- 
toevaporation by UVB is not complete and the effects of 
self-shielding are still significant. Hence, in the context 
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of CDM cosmology, some fraction of baryons forms stars 
inevitably. Thus, we have to take into account the shelf- 
shielding properly, in order to evaluate the effects of UVB 
feedback on substructure problem. 

To make a further approach on this problem, other ef- 
fects such as the internal UV sources, the evaporation 
driven by SN explosions with the strongly top-heavy ini- 
tial mass function of formed stars or the gas tripping in 
intracluster medium may have to be investigated. 

4.2. dE/dSph's in Local Group 

What objects correspond to the simulated galaxies? 
One possibility is dE/dSph galaxies in the Local Group 
(LG). The dwarf galaxies in LG are known to be extended 
(low surface brightness), spheroidal system (e.g. Mateo 
1998). Especially, faint dSph's (My > -14) have a larger 
mass-to-light ratio. Hirashita, Takeuchi, & Tamura (1998) 
have found 

log(Mvir/i) = 2.01ogMvir-hl3.1 (16) 

for Mvir < IO^Mq, while log(Mvir/i) ^ 0.4 - 0.7 for 
Mvir > W^Mq. These observational features indicate that 
the mass loss is more prominent for Mvir < lO^M©. This 
can be attributed to SN-driven winds (Hirashita, Takeuchi, 
& Tamura 1998; Mori, Ferrara, & Madau 2002). The 
present numerical simulations suggest that such a wind 
can occur also by the photoheating due to the reioniza- 
tion, even before SN explosions drive a thermal wind. If 
wc use the final stellar fraction (14), we find Mhaio/^ = 63 
with assuming M^tai/L = 3 for Zc = 1-7 and Mhaio — 
6 X IO^Mq. Such a high M/L is also obtained by the re- 
cent results by Ricotti (2003) in which 3D cosmolgical 
radiation hydrodynamical simulations arc performed with 
OTVET approximation (Gnedin & Abel 2001). We also 
remark that Shaviv & Dekel (2003) pointed out the con- 
nection between dSphs and photoevapolated galaxies with 
ID simulations. 

In addition, the SFH in the Local Group dwarfs has a 
ubiquitous feature (e.g. Mateo 1998; Gnedin 2000a). The 
stars in most LG dwarfs are estimated to have formed 
more than 10 Gyrs ago. This means that star formation 
rate should have a clear cut-off at such epoch. If we as- 
sume the age of the universe is 13 Gyr, 10 Gyr in look- 
back time corresponds to 2: ^ 3, which is 2 Gyrs later than 
the cut-off epoch (2: ~ 4 — 5) in the present simulations. 
However, the observed star formation history itself should 
have uncertainties of a few Gyrs. Thus the reheating by 
reionization could be one of the promising mechanism to 
provide such a cut-off in SFH of the LG dwarf galaxies. 
Also, the SFH in the LG dwarfs has another important 
feature: there still remain weak star formation activities 
after the cut-off of star formation rate. This feature may 
be reproduced by the numerical galaxies presented here, if 
c* is at a level of 0.01 for relatively high mass and high-2;c 
models (see §3.2). Furthermore, if one incorporates the re- 
cycling of gas ejected from the stars, which is not included 
in the present simulation, additional stars may form unless 
the compete evaporation does not occur. This tendency 
is expected to be conspicuous for the halos with a deep 
gravitational potential. However, the gas ejected to the 
intersteller space in such dwarf galaxeis will be photoion- 
ized again. The effects of self-shielding and the formation 
of dust will delay the photoevaporation of the processed 



gas and consequently some amount of gas will be converted 
into stars again. 

4.3. dirr's in Local Group 

On the other hand, the observed dirr's in LG still have 
significant star formation activities embedded in old ex- 
tended spheroidal components. Although we have no 
counterpart in the present simulations, we can infer the 
origin of such galaxies by extrapolating the numerical re- 
sults. Massive (~ lO^^M©) halos formed at low Zc might 
be the candidate of such galaxies. In such halos, small 
scale perturbations grows slowly. Thus, the fraction of 
gas component after the reionization of the universe is 
largo as seen in Figure 7. In addition, suc;li halos have 
a gravitational potential deep enough to sustain photoion- 
ized gas, so that stars form efficiently from such abundant 
gas. Therefore, active star formation is expected even af- 
ter the reionization. We can suggest that the conditions 
requisite for dirr formation are 1) that the gravitational 
potentials of dark halos are deep enough to retain photo- 
heated gas, and 2) that the galaxy formation epoch is later 
than the reionization epoch. Actually, dirr's appear to be 
rather massive (e.g. Mateo 1998). Also, the reduction of 
star formation rate due to UVB may suppresses the SN- 
driven galactic wind and thereby allows the formation of 
diffuse dwarf systems like Irr's. 

5. CONCLUSIONS 

In this paper, we have performed 3D radiation hydro- 
dynamic simulations on the formation of low-mass objects 
under UVB. The suppression by UVB of the formation of 
low-mass objects at lower redshifts (M < lO^M©, Zc < 5) 
is confirmed. But, simultaneously it is found that the for- 
mation of low-mass galaxies at low redshifts is not com- 
pletely forbidden by the UVB-induced photoevaporation. 
Baryons at high-density peaks can collapse and be self- 
shielded from UVB even during the reionization, to form 
stars eventually. This occurs even if the universe was reion- 
ized at earlier epochs as suggested by the WMAP. Thus, 
effects of UVB feedback on the substructure problem with 
a cold dark matter (CDM) scenario should be evaluated 
with careful treatment of the radiative transfer. 

In order to clarify the effects of radiative transfer, we 
have compared the results of simulations with those under 
the assumption of optically-thin medium. 

The reduction for the final stellar fraction is 30-40% 
without radiative transfer in the ^Jreion — 7 regime. Hence, 
the star formation before the reionization is important to 
determine the final stellar fraction. But, the radiative 
transfer effect is much more serious in the earlier reioniza- 
tion (^reion — 17) regime, where the reduction for the final 
stellar fraction is more than an order of magnitude. Also, 
it should be stressed that in the optically-thin approxima- 
tion, gas temperature of the cloud becomes always higher 
than 10** K after the reionization, so that stars cannot 
form there, whereas the self-shielding arising from radia- 
tive transfer allows star formation even after the reioniza- 
tion if the gravitational potential is deep enough to retain 
the gas. The star formation after the reionization must 
influence the subsequent evolution of dwarf galaxies with 
the possible recycling of interstellar matter. 
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In the previous optically-thin simulations (e.g. Tassis 
et al. 2003), the star formation criterions are different 
from those in the present simulation. The conditions used 
here are stricter than those used in the previous works. In 
previous simulations, stars can form slowly even after the 
reionization due to the milder star formation conditions, 
although the cold HI regions cannot be formed. Although 
such a treatment of star formation might be approximately 
valid, it should be confirmed by the simulations with ra- 
diative transfer. 

We also have investigated the effects of c*, and found 
that the final stellar fractions are not so different for var- 
ious values of c, . However, a smaller value of c, delays 
star formation and allows enhanced star formation after 
the reionization if the system is relatively massive and the 
collapse epoch is considerably later the reionization. Also, 
if the stellar feedback is taken into account, the further 
effects of c* may be expected, since 1) larger results in 
faster recycling (i.e. a larger amount of heavy elements), 
and 2) larger c* is followed by simultaneous SN explosion, 
which means easier disruption of host galaxies. These is- 
sue should be investigated by the simulations that take 
into account the stellar feedback coupled with radiative 
transfer. 

The star clusters formed at high-density peaks coalesce 
with each other in a dissipationless fashion to form a 
spheroidal system with a large mass-to-light ratio. It is 
also found that there is a sharp cut-off of star forma- 
tion rate at z ~ 4 — 5 in all the simulations because 
of the appreciable photoevaporation of cold gas clouds 
in the reionized universe. Moreover, in massive dwarf 
galaxies (Mhaio — 6 x lO^M©) formed at high redshifts 
(zc ^ 8), weak star formation ac;tivities continue even 
after the reionization. On the other hand, less massive 
(.^haio ^ 6 X IO'^Mq) galaxies formed later {zc < 8) do 
not have star formation activities after z ~ 4 — 5. Obser- 
vational counter parts of these systems might be the dwarf 
spheroidal galaxies in the Local Group. Those galaxies 
have extended, old and spheroidal stellar systems, with 
characteristic star formation histories. It is also known 
that fainter galaxies have larger mass-to- light ratio. All 
of those features agree qualitatively with those of the sim- 
ulated galaxies, although the stellar feedback should be 



also investigated to make quantitative comparison with 
observed galaxies. 
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APPENDIX 
TESTS FOR NUMERICAL SCHEME 

Here, we show some tests of the present numerical scheme. 
Using the above scheme, we tested the propagation of ion- 
ization front in a uniform media. The test results are 
shown in Figs. 13 and 14, in which several snapshot of 
the ionization structures and temperature distributions are 
shown. In this test, uniform cube is irradiated from the 
left by the ionization flux which is proportional to u^^. 
The initial (neutral) optical depth of each SPH particle at 
the Lyman limit is ~ 10^. The vertical axis denote the 
fraction of electrons (Fig. 13) / temperature (Fig. 14), and 
the horizontal axis shows the normalized size of the cube. 
The position of the particles are located at the grid points. 
The axes of the grid are set so that they are not parallel 
or perpendicular to the surfaces of the cube. The vertices 
denote the ionization structure at four different times by 
the RSPH code, and the four solid lines denote the results 
obtained by ID mesh code, which uses 2 x 10''' meshes (this 
number of mesh corresponds to (5t ~ 0.1 at Lyman limit, 
that is enough to resolve the ionization front). It is quite 
clear that RSPH scheme can reproduce the propagation of 
ionization structure fairly well, although the optical depth 
of an SPH particle is much larger than unity. 
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Table 1 



Properties of Formed Galaxies 





Mhalo/M0 


Tvir/K 


CTiD/km s ^ 


criD(DM)/km s-i 


-Rhalf/pC 


/star 




1.7 


5.9 X 10^ 


3.6 X 10^ 


4.6 


5.5 


62 


0.28 


63 


3.4 


5.9 X 10^ 


5.1 X 10^ 


5.1 


6.5 


58 


0.37 


50 


5.1 


6.1 X 10^ 


6.3 X 10^ 


5.7 


7.2 


69 


0.47 


44 


6.5 


6.1 X 10^ 


8.2 X 10^ 


6.4 


8.2 


66 


0.61 


36 


8.1 


6.1 X 10^ 


1.0 X 10"^ 


7.1 


9.2 


62 


0.71 


33 


1.1 


6.0 X 10* 


1.7 X 10"^ 


14 


12 


2.9 X 10^ 


0.43 


40 


2.9 


6.0 X 10* 


2.0 X 10"^ 


16 


14 


2.1 X 10^ 


0.55 


31 


4.5 


6.0 X 10* 


2.0 X 10^ 


17 


16 


2.2 X 10^ 


0.68 


31 


6.0 


6.1 X 10* 


2.1 X 10"^ 


19 


19 


1.8 X 10^ 


0.79 


29 


7.6 


6.1 X 10* 


2.6 X 10* 


20 


21 


1.6 X 10^ 


0.93 


26 
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Fig. 1. — H2 Cooling functions induced by H-H2 collision are plotted against the temperature for two different densities. Solid curves axe 
plotted based upon the formula given by HoUenbach & McKee (1979), long dashed curves by Galli & Palla (1998) and the short dashed 
curves by Lepp & ShuU (1983). They are discrepant at low density. In this paper, those by Galli & Palla (1998) are adopted. 




Fig. 2. Cooling functions below lO^K arc plotted. Two solid lines denote the cooling rate due to metals, assuming Z/Zq = 1 and 10~^ 
respectively. H2 cooling function is plotted for the postshock region with shock velocity larger than 30 km/s {dotted line) and also for a given 
abundance of 2/H2 = 10~* {short dashed line). A long dashed line represents the rate by hydrogen Ly a. 
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(1 ) Projection onto x-y plane (2) Domain decomposition on x-y plane 




Z 

Fig. 3. — Parallclization of the code is schematically shown. (1) Particles arc projected onto the x — y plane which is perpendicular to the 
light ray. (2) The particles are sorted by x— and y— coordinates and decomposed into subclasses which include an equal number of particles. 
(3) Each subclass is assigned to each PU. 
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Fig. 4. — Time evolution of the spatial distributions of SPH and star particles for A/ijaio — 6 X IO'^A/q and Zc ~ 1.7 is shown. The figure 
consists of 8 panels. Upper-left panel is the distribution of SPH particles at very early phase, and it evolves downward in the left column, 
continuing to the upper-right panel. Corresponding redshift is printed at the upper-left corner of each panel. The colors of dots represent the 
logarithmic temperature of the SPH particles, and the color legend is shown at the bottom in logarithmic scale, as well as the physical scale 
ruler (3 kpc) at the right-bottom corner. The blue particles do not represent SPH particles, but star particles. 
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Fig. 5. — Same as Fig. 4 except that the parameters are A/haio — 6 X 10* Af© and Zc ~ 7.6 
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Fig. 6. — The time-variation of star formation rate (SFR) is plotted against the cosmic time. Four panels represent respectively the results 
of four different runs, where the parameters are displayed in each panel. The vertical short dashed lines denote z = 7 and z = 6, between 
which the reionization is almost completed in this model. 




Fig. 7. — The final stellar fraction /star is plotted against the collapse rcdsliift. Each mark denotes the result of a run: pentagons(c, = 1), 
squaj:es(c« = 0.1), and triangles(c* = 0.01) denote the runs with Mhaio — 6 x IO'^Mq. Similarly, starred pentagons(ct = 1), starred 
squaj:es(c* = 0.1), and starred triangles(c* = 0.01) represent the results of relatively massive galaxies (Mhaio — 6 x 10 Mq). 
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Fig. 8. — Same as Fig. 6, except that c» = 0.01. 




Fig. 9. — The time-variations of stellar fraction /star arc plotted against the cosmic time for transfer and optically-thin cases. Solid lines 
are the results with radiative transfer, while dashed lines arc those under the assumption of optically-thin. Parameters of four panels are the 
same as that in Fig.6. The vertical short dashed lines denote z = 7 and 2 = 6, between which the reionization is almost completed in this 
model. 
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Fig. 10. — The final stellar fraction /star is compared between radiative transfer and optically-thin cases, assuming c* = 0.1. Filled symbols 
represent the results of full radiative transfer simulations, and open symbols denote the results of optically-thin simulations. 
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Fig. 11. — The rotation energy-to-random motion energy ratio (i?rot/i?ran) is plotted against the cosmic time. Parameters of four panels 
are the same as that in Fig. 6. Spikes in this figure come from merger of subclumps. Since Brot and Bran are measured from the center of 
gravity (see the text), violent mergers cause high spikes. Especially, in the left panels, the collapse epoch is relatively late, so that such violent 
merger phases are conspicuous. 
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Fig. 12. — Star formation histories in the early reionization regime are shown (upper two panels), and time evolutions of stellar fraction 
are shown by the solid curves (lower two panels). The parameters are attached in each panel. The results of optically-thin simulations are 
also shown in the lower two panels (dashed curve). In these simulations, = 0.1 is assumed. The vertical short dashed lines denote z = 17 
and z=16, between which the reionization is almost completed in this early reionization model. 




Fig. 13. — Propagation of ionization front is shown (see also the text). Crosses denote the numerical results, while solid lines are analytic 
solution. 
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